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We study the variation of tlie dielectric response of a dielectric liquid (e.g. water) when a salt is 
added to the solution. Employing field-theoretical methods we expand the Gibbs free-energy to first 
order in a loop expansion and calculate self-consistently the dielectric constant. We predict analyt- 
ically the dielectric decrement which depends on the ionic strength in a complex way. Furthermore, 
a qualitative description of the hydration shell is found and is characterized by a single length scale. 
Our prediction fits rather well a large range of concentrations for different salts using only one fit 
parameter related to the size of ions and dipoles. 



Electrostatic interactions in aqueous media between 
dipoles and charged objects such as ions, colloidal par- 
ticles and interfaces, play an important role in electro- 
chemistry, biology and materials science. The century- 
old Poisson-Boltzmann (PB) theory gives a simple and 
powerful description for numerous systems, taking into 
account only Coulombic interactions on a mean-field 
level, treating ions as point-like particles and the aque- 
ous solution as a continuous and homogeneous dielectric 
medium with a constant dielectric constant, ew The PB 
theory succeeded over the years in capturing much of the 
underlying physics of electrolyte solutions, and, in partic- 
ular, is successful when applied to monovalent ions and 
weak surface charges pJ-lJl- 

The PB theory has several limitations. It does not take 
into account neither the correlations between the charges, 
nor does it allow for any fluctuations beyond mean field. 
This leads to significant unaccounted corrections in cases 
of high charge density, especially near charged surfaces 
and interfaces Q . In order to improve upon the PB the- 
ory, several extensions have been offered in recent years, 
and effects of correlations and fluctuations for charged 
interfaces were considered, for example, by using inte- 
gral equation theories @. In other approaches 0], the 
PB theory was modified in a simple and elegant way to 
include steric and other ionic-specific effects preventing 
ions from accumulating near charged surface for very high 
ionic concentrations. Furthermore, in the well known 
Deryagin-Landau-Verwey-Overbeek (DLVO) theory 0, 
van der Waals attractive interactions were added to the 
electrostatic repulsion in order to explain charged col- 
loidal stability. More recently, Monte-Carlo (MC) sim- 
ulations and Molecular Dynamics (MD) were employed 
in order to study specific solvents and solutes and the 

interaction between them urn. 

The PB theory as well as others primitive models of 
ionic solutions |1J] assumes that the ions are immersed 
in a continuum dielectric background characterized by 



the dielectric constant of water, e^,. Hence, in order 
to model the experimentally known dielectric decrement 
phenomenon [3, [l5l - [l7l | , one needs to employ more refined 
theories, which take into account ion-dipole correlations 
and fluctuations. The electric decrement stems from the 
fact that the local electric field around each ion is greater 
than the external field, and orients the dipolar water 
molecules in its vicinity. This creates a hydration shell 
[lit i^Skt encircles the ions. The total response of dipoles 
to the external field is thus smaller and leads to a re- 
duction in the dielectric constant. In dilute solutions the 
dielectric constant depends linearly on c,,, the salt con- 
centration, e{cs) = Ew — "fCs, where 7 is the coefficient of 
the linear term. The value of 7 is ion specific and ranges 
from 8M"i to 20M"\ for concentrations up to 1.5 M. 
At higher c, values, noticeable deviations from linearity 
are observed and are due to ion-ion interactions [ISE 

In this Letter we go beyond the standard PB theory 
and investigate the microscopic origin of the dielectric 
decrement. We construct the field theory for a grand- 
canonical ensemble of point-like ions and dipoles. On 
a mean-field level, we obtain an analytic description of 
the hydration shell around ions and calculate the aver- 
age dielectric constant in the shell. Then, using a field- 
theoretical approach we examine the dielectric decrement 
by a one-loop calculation. A closed-form formula for the 
dielectric constant is obtained and depends on a single 
length-scale, related to the typical size of ions and water 
molecules. Our results show a substantial contribution 
arising from the one-loop correction terms, and deviates 
from the linear relation, e(cs) = Em — jCs- We fit the 
calculated e(cs) and compare it with the experimental 
data for several monovalent ions. The results fit rather 
well for a wide range of ionic concentrations, even in the 
non-linear decrement regime. 

Although our formulation is very general, in this Let- 
ter we focus on the dielectric constant variation for ionic 
solutions. The solvent (water) is modeled as a liquid of 
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point-like permanent dipoles of dipolar moment po and 
the symmetric 1:1 monovalent salt as a liquid of point-like 
charges, ±e. Using the standard Hubbard- Stratonovich 
transformation [l9[, the grand-canonical partition func- 
tion Z and free energy F = logZ[p{r)] in presence 
of an external charge distribution p{r) can be written as 
a functional integral over the field 0(r) that is conjugate 
to p(r) 



Z[p(r)] ^ / l?0(r)cxp -/3 / d^r /((/.(r)) 



t(3 / d3r(/)(r)p(r; 



(1) 



with 



A, 



[V</.(r)]2- 2A, cos[/?e0(r) 
sin(/3po|V^(r)|) 



(2) 



/3po|V0(r)| 

where 13 = l/ksT, and /cbT is the thermal energy. The 
fugacities of the salt and water A^ are determined by 
solving the implicit equations 

V dXs 

Xd_d_ 

V dx. 



■ log Z = , 

■ log Z ^ cd , 



(3) 



where V is the volume and Cg and are, respectively, 
the salt concentration and density of water molecules. 
The electrostatic potential ^'(r) is given by 



vl/(r) 



1 ^logZ[p(r)] 



(4) 



where the bracket (...)o denotes the thermal average with 
the Boltzmann weight given in Eq. ([1]), for p = 0. 

The Gibbs free energy G is a function of ^'(r) and is 
defined as the Legendre transform of F[p{r)] with respect 
to p{r) through the equations 



G[^{r)] = F[p{v)] 
6F[p{r)] 



d'V *(r)p(r) 



5p{r) 



from which we deduce the Legendre relation 



P(r) 



<S*(r) 



(5) 



(6) 



The dielectric tensor Ea/B is defined through the Fourier 
transform: 



dpadpi3 



(52G[*(r)] 



(7) 



*=0,p=0 



In our present study, the dielectric tensor for isotropic 
aqueous solutions is diagonal e^^ = eSap, and our aim is 



to calculate its variation as a function of the salt concen- 
tration Cs- Following Eq. ([7]), we need only to calculate 
the coefficient of the |V\I'(r)p terms in the Gibbs free- 
energy G[4'(r)]. 

On the mean-field level the Gibbs free-energy for the 
Dipolar PB system, Gdpb, is 



/3Gdpb[*] 



' Cd9{u), 



2cs cosh(/3e*) 



(8) 



with u = /3poV5'(r), u = |u|, and g{u) = sinhu/u. The 
Dipolar Poisson-Boltzmann (DPB) equation is an exten- 
sion of the standard PB equation, and can be derived via 
a variational principle from Eq. ([5]): 



eoV^* = -2csesinh {l3e^) 



e(/3po|Vvl/| 



(9) 



where the function Q = g'{u) ~ coshu/u — smliu/v? is 
related to the Langevin function L{u) = coth(u) — 1/u. 

The dielectric constant e^PB is calculated by substi- 
tuting GppB into Eq. ([7]). The result is the same as in 
Ref. 13: 



(10) 



As can be seen above, the dielectric response depends on 
the dipole (but not the ion) density Cd and, hence, cannot 
explain any dielectric decrement. At room temperature, 
T = 300 K, and for pure water with dipolar moment pq — 
1.8 D, and density Cd = 55 M, the obtained value of e is 
e-^pg ~ 11.1. Note that this value is much smaller than 
the measured water value Sw — 80. This is not surprising 
since Eq. (jlOp is a dilute gas approximation and does not 
capture the correlation effects in the case of dense liquid 
water as well as the finite-size of water molecules. 

It is still possible, however, to circumvent this draw- 
back by solving the DPB equation, Eq. around a 
fixed point-like ion, and showing the existence of a hy- 
dration shell around it. Expanding ^ in a Taylor scries, 
the dielectric constant is extracted as a function of the 
distance r from the ion: 



e(r) 



3h^ ilhlr) + 1 ' 



(11) 



where the function h{u) is obtained by solving a cubic 
equation: 



h{u) 



1 u -/u 
27 ^ 4 ~2~ 



1/3 



1 U y/u 

27 ^ 4 ^ T 



-,1/3 



(12) 



The length, ll = lBd/^/\Q in Eq. depends both 

on the dipole size d — pa/l and the Bjerrum length. 
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Ib = / A-Kek-oT . It is the only length scale in the prob- 
lem characterizing the hydration shell thickness as can 
be seen from the behavior of e(r) (inset of Fig. 1). In 
the vicinity of the ion (r < Ih)^ the dielectric response 
is very small and it smoothly rises to bulk values as the 
influence of the ion decreases, within distances of a few 
Ih- 

By averaging e of Eq. ([TT|) over a sphere of radius R 
around each ion, and equating R with the typical dis- 
tance between two ions, (2cs)~^/^/2, at concentration 
Cs, we obtain the expression of {e{cs)). As can be seen in 
Fig. 1, the non- linear dielectric decrement is reproduced 
and fits rather well the experimental data for RbCl and 
CsCl salts of Ref. For low salt concentration, the 

averaging of e can be done analytically and results in 
7 ~ llOe^pgZ^Cs, which is proportional to the hydration 
shell volume. With parameter values as in Fig. 1, we get 
7 ~ 17. 




FIG. 1. (color online). The dielectric constant (e) averaged 
inside a specific volume around a single ion (solid line) as func- 
tion of ionic concentration, Cs. The comparison is done with 
experimental values for RbCI (empty circles) and CsCI (full cir- 
cles) [16]. In the inset, the exact (solid line) and approximated 
(dashed, Eq. ([11])) solutions of the DPB equation, Eq. dO} are 
shown as function of the distance r from a point charge (ion). 
Choosing as a fit parameter the dipole moment of water to be 
po = 4.8 D (instead of the physical value po = 1.8 D) [2Q|]. allows 
us to obtain e^pB = 80 and lu — 1.5 A. 

So far in order to induce the necessary ion-dipole cor- 
relations, we had to rely on calculating the dielectric re- 
sponse around a single fixed ion, which gives the dielec- 
tric decrement in an approximated way. Moreover, the 
validity of the DPB theory can be justified only in the 
dilute limit, where the dipole-dipole fluctuations are not 
important. 

To overcome this limitation a more complete treat- 
ment of the statistical mechanics of ionic and dipolar 
degrees of freedom (including their interactions) is de- 
veloped. We generalize the Debye-Hiickel approximation 
to dipolar systems in presence of salt, via a loop expan- 



sion [211 l22l of the Gibbs free-energy, G = Gdpb + AG. 
To one-loop order the calculation of AG yields: 

AG[*(r)] ^Trlog [-eoV^ + 2A,s/3e2 cosh(/3e*(r)) 
+ \dPpl {dj:^^{v)d^ + F„^(r)9„ • d^)] (13) 



with 



r.,(r) = 5„,^-!^f^-,» 



(14) 

To get a consistent expression for the dielectric con- 
stant we need to pay extra attention in the loop expan- 
sion to the expressions of the ion and dipole fugacitics. 
As and A^. While at mean- field level As and Ad are equal 
to the ion and dipole (water) densities, Cs and c^, respec- 
tively, their corrections are derived by expanding Eq. ([3]) 
to one-loop order. To this order, only the correction to Ad 
will affect the mean-field value of the dielectric constant. 



ej3P3 , and is written as 



Cd 



2tt Sr- 



--0 



3a3 



1 - 

3 , ,Q _i 27r 
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where a is a microscopic cutoff length and kd = VStt^bCs 
is the inverse Debye length. The origin of the cutoff in 
our formulation is related to the fact that Coulombic in- 
teractions diverge at zero distance, while in reality such 
a divergence is avoided because of steric repulsion. A 
simple way to take this into account is through a cutoff 
length a, which is related to the minimal distance be- 
tween adjacent dipoles and charges, and thus indirectly 
also to the size of dipoles and ions. 

We can now calculate consistently the corrections to 
the dielectric constant up to one-loop order. The correc- 
tion term can be split into water and salt contributions. 



Aed = 
A£s = 



(^DPB ^ 



47r 



(^DP 



3cda3 



1 tan 

27r 



271 



(16) 



The term Asd represents the fluctuation effect of the wa- 
ter dipoles beyond the mean-field DPB level. It varies as 
^ l/(cda^). This pure water fiuctuation term essentially 
adds a positive numerical prefactor of rather large mag- 
nitude to the mean- field value of e^PB ' nieaning that the 
one-loop correction is important even for the pure water 
case. 

The second term Aes has by itself two contributions. 
The leading term in the dilute solution limit, kdo ^ 1, 
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depends linearly on the salt concentration, Ae^ = —70^ 
with 



^DPB 



(17) 



When the Debye length ^ is of the same order of mag- 
nitude as a (high salt limit), the last term in Ae^ starts to 
dominate and the dielectric decrement becomes smaller 
until eventually it will reverse the trend and may even 
cause a dielectric increment as seen in experiments [l6j 
in the high salt limit. 

We compare our prediction for the dielectric constant 
£, Eq. p^ . to experimental values [31 for seven different 
ionic solutions in concentration range of 0-6 M. We sepa- 
rate the seven salts into three subgroup according to the 
size of the alkaline cations as presented in the three parts 
of Fig. 2. In each of the figure parts the a parameter is 
fitted separately. We treat a as a free parameter and find 
its value by the best fit of our prediction, Eq. (fT6|) , to ex- 
perimental data, while keeping the physical value of the 
water dipolar moment, po = 1.8 D. The best fit to the 
data can be seen in 2(a) and corresponds to the largest 
ionic size of Cs"*" and Rb+. In (b) the fit for K+ ions is 
also quite good, while in 2(c) for the smallest ions, Na"*" 
and Li"*", the deviation is more pronounced, especially for 
the larger values of Cg > 4M. 

Note that our formula takes into account only in a 
broad sense the finite size of ions (and the distance of 
closest approach between them) via the a parameter. It 
is beyond the level of the theory to give more specific 
ionic predictions. Hence, the obtained value of a ~ 2.7 A 
is not very sensitive to the type of salt, but its main con- 
tribution comes from the water dipoles themselves whose 
diameter is about 2.75 A. On the other hand, as can 
be clearly seen from Fig. 2, important cooperative effects 
of ions and dipoles are accounted for in our non-linear ex- 
pression for e{cs). For small Cg, the dashed line represents 
the best linear fit and works well only when Cs < 1 M, 
while the non-linear prediction of Eq. (jl6p succeeds in 
fitting the large concentration range as well [23 |. 

In conclusion, we are able to reproduce rather well the 
linear and non-linear dielectric decrement behavior over a 
large range of ionic concentrations, and the obtained val- 
ues of e are in quantitative agreement with the data for 
several types of monovalent salts. In addition, we found 
a qualitative description of the hydration shell charac- 
terized by a single length scale, Ih- Note that our model 
does not contain any significant ionic-specific effects. To 
improve on the latter, we would need to include the ionic 
finite-size and specific non-electrostatic short-range inter- 
actions that affect both bulk properties of ionic solutions 
such as dielectric constant and viscosity as well as their 
behavior at interfaces and, in particular, their surface 
tension. 
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FIG. 2. (color online). Comparison of the predicted dielectric constant, e, from Eq. (|16p with experimental data from Ref [l^ . 
as function of ionic concentration, Cs, for various salts. The theoretical prediction (solid line) was calculated using a as a fitting 
parameter. In (a) the fit is for RbCI and CsCI salts with a = 2.66 A; in (b) the fit is for KF and KCI with a — 2.64 A; while in (c) the 
fit is for LiCI, Nal and NaCI with a = 2.71 A. The dashed lines are the linear fit to the data in the low Cs < 1 M range. The slope of 
the linear fit is 7 = 11.7 M"^ in (a), 9.0 M"^ in (b) and 13.8 M"^ in (c). The value of 7 for each salt varies by about 10-20 % and 
the linear fit should be taken as representative for the low Cs behavior. 



